TECHNICAL  REPORT  NO.  577 


INVISCID  FLOW  FIELD  PAST  A  POINTED 


CONE  AT  AN  ANGLE  OF  ATTACK 


PART  I  -  ANALYSIS 


By  G.  Moretti 


C  L  E  A  R  I  N  G  li  0  U  3  t 
POR  HT.DBRAL  SC’IKN^ !  IFIC  AND 
TECHNICAL  1NI‘T>KMATI0N 
HBnaoopxn^  Miorofithal  j 


^  } 

^  / 1 


I  y.'  Y/ 


:>  PP! 


Cf-rC.  / 


December  15,  1965 


■Mfiilv 


Project  8030/0650 


Total  No.  of  Pages  -  iv  and  40 
Copy  (12)  of  125 


TECHNICAL  REPORT  NO.  577 

INVISCID  FLOW  FIELD  PAST  A  POINTED 
CONE  AT  AN  ANGLE  OF  ATTACK^ 
PART  I  -  ANALYSIS 

By  G -  Moretti 


Prepared  for 

Advanced  Research  Projects  Agency 
Washington,  D.  C. 

Under  Contract  SD-149,  Supplement  No.  6 
ARPA  Order  No.  396 
Project  Code  3790 

"Ballistic  Reentry  Studies" 

Project  Engineer  -  Dr.  H.  Lien 
Code  516  -  333-6960 


Prepared  by 

General  Applied  Science  Laboratories,  Inc. 
Merrick  and  Stewart  Avenues 
Westbury,  L-  I.,  New  York 


December  15,  1965 


Approved  by:- 


*This  research  is  sponsored  by  the 
Advanced  Research  Projects  Agency. 


Robert  W.  ByrAle 


Vice  Presidi^ 


tiiMpp- 


XI 


SUMMARY 


The  conical  inviscid  flow  field  past  a  circular  cone  at  an 
angle  of  attack  is  obtained  as  the  asymptotic  stage  of  a  three- 
dimensional  flow.  The  latter  is  computed  in  a  spherical  frame 
by  means  of  finite  difference  techniques  associated  with  a  method 
of  characteristics  on  the  shock  and  body.  A  brief  critical  dis¬ 
cussion  of  previous  papers  on  the  subject  and  an  analysis  of  the 
numerical  results  are  given.  Comparison  with  experiments  shows 
that  the  applicability  of  the  theory  is  severely  limited  by 
boundary  layer  effects. 
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TECHNICAL  REPORT  NO.  577 
INVISCID  FLOW  FIELD  PAST  A  POINTED 
CONE  AT  AN  ANGLE  OF  ATTACK 
PART  I  -  ANALYSIS 

By  G.  Moretti 

I.  INTRODUCTION 

The  State  of  the  Art 

The  steady  supersonic  flow  past  a  pointed  cone  of  revolu¬ 
tion  at  zero  angle  of  attack  was  first  studied  by  Busemann  in 
1929,  and  by  other  authors  in  the  thirties  (Ref.  1).  The 
definition  of  a  conical  flow  was  given  as  that  of  a  flow  whose 
parameters  do  not  depend  on  r  in  a  spherical  frame  of  refer¬ 
ence,  (r,0,(p).  As  a  matter  of  fact,  the  flow  past  a  pointed 
cone  of  revolution  at  zero  angle  of  attack  is  also  axisym- 
metric,  that  is,  its  parameters  do  not  depend  on  <p  either.  The 
latter  problem  is  reduced  to  solving  two  ordinary  differential 
equations,  with  0  as  the  only  independent  variable; 

f  du/d0  -  V 


dv/d0  == 


u  +  (u  +  V  cot  0)/(v“/a~  -  1) 


(1) 


2 

(the  three  velocity  components  in  the  spherical  frame  being 
U/  V,  and  w,  respectively,  and  a  being  the  speed  of  sound) . 

In  1947,  such  equations  were  coded  for  an  electronic  computer 
by  Kipal  and  his  staff  (Ref.  2) ,  and  the  results  tabulated. 
Today,  Kopal's  tables  are  obsolete,  since  a  program  of  not 
more  than  20  FORTRAN  statements  may  provide  a  direct  evalu¬ 
ation  of  the  flow,  once  the  free- stream  Mach  number,  the 
semi-angle  of  the  cone,  e,  and  the  ratio  of  specific  heats, 
y,  are  given. 

In  the  late  forties,  the  problem  of  a  pointed  cone  at 
an  angle  of  attack,  Ot,  was  considered,  with  some  optimism, 

t 

to  be  a  simple  extension  of  the  case  above,  at  least  for 
®  It  is  easy  to  prove  that  a  steady  conical  flow  may 

exist,  whose  parameters  depend  on  0  and  (p.  .  Five  differ¬ 
ential  equations  are  needed  now  (since,  in  addition  to  u 
and  V,  the  third  velocity  component,  w,  the  pressure,  p,  and 
either  the  density,  p,  or  the  entropy,  S,  must  be  detemined)  , 
but  the  independent  variables  number  only  two.  In  fifteen 
years,  several  attempts  have  been  made  to  solve  the  problem 
so  formulated,  and  the  number  of  papers  is  so  la  rge  that  an 
exhaustive  list  cannot  be  given  here. 


ih. 
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The  problem  proved  to  be  a  very  tough  one,  since  the 
system  of  equations  which  describes  it  is  elliptic,  if  ot 
is  small  and  v®+w*<a*,  and  becomes  of  the  mixed  type  for 
larger  values  of  Ot,  if  there  are  regions  of  the  ( 6, fp) -sphere 
in  which  v*+w^a®.  In  such  a  sphere,  the  shape  of  the  cone 
is  given,  but  the  shock  wave  is  unknown.  In  this  respect, 
the  problem  of  a  cone  at  an  angle  of  attack  is  similar  to 
the  blunt  body  problem,  and  similar  techniques  have  been 
used  in  both  problems,  such  as  inverse  methods  with  step-on 
integration  from  the  shock  to  the  body  (Ref.  3),  methods 
of  integral  relations  (the  "strip-method"  of  Belotserkovsky) 
(Ref.  4),  and  methods  of  series  expansions  (Ref.  8).  In  the 
latter  case,  some  papers  try  to  improve  the  original  linear¬ 
ized  formulation  of  Stone  (Refs.  5,6)  and  Ferri  (Ref.  7), 
others  (Refs.  9,10)  are  related  to  Chernyi's  expansion  in 
powers  of  (y-l)/(y+l),  etc. 


~ifcr 
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II.  THE  CENTRAL  IDEA  FOR  A  NEW  TECHNIQUE 

Critical  objections,  similar  to  the  ones  discussed  in 
some  detail  for  the  blunt  body  problem  (Ref.  11) ,  may  be 
raised  here.  Briefly,  since  a  numerical  procedure  is  used 
to  solve  the  problem,  a  system  of  elliptic  equations  is 
ill-suited  for  the  job.  In  the  blunt  body  problem,  a 
powerful  and  accurate  technique  may  be  set  up  if  the  system 
is  made  hyperbolic  by  considering  the  steady  state  as  the 
asymptote  of  an  unsteady  one.  In  the  present  problem,  the 
equations  of  a  steady  flow  in  the  three  coordinates,  r,  6, 
and  <p,  fom  a  hyperbolic  system,  consequently,  the  problem 
is  easier  to  handle  if  the  motion  is  not  conical.  Let  us 
consider,  then,  the  conical  motion  as  the  asymptote  (with 
respect  to  r)  of  a  non-conical  one.  To  this  effect,  assume 
an  arbitrary  set  of  initial  conditions  on  a  sphere,  r  =  r^, 
say,  and  compute  the  flow  from  there  on  as  dependent  on  0, 

<p,  and  r,  until  an  asymptotic  state,  independent  of  r, 
develops.  The  physical  argument  for  the  existence  of  a 
conical  flow  as  an  asymptotic  one  is  that  the  flow  in  the 
shock  layer  is  forced  by  the  outer  flow,  which  is  uniform 
(and  therefore  conical)  and  the  geometry  of  the  body  (which 
is  a  cone) . 


tU- 
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In  the  present  approach  an  initial  value  of  r  is  defined, 
and  we  suppose  that  the  asymptotic  values  are  reached  at 
another  value  of  r,  r  =  ri ,  say.  The  difference,  rj  -  r^, 
represents  the  actual  length  which  the  flow  should  travel  in 
order  to  overcome  the  initial  state  (a  perturbed,  non-conical 
one)  with  a  continuous  build-up  of  conical  effects.  Such  a 
length  will  obviously  be  the  longer,  the  farther  the  initial 
conditions  are  from  a  conical  flow.  But  the  length  itself  is 
immaterial,  since  no  length  is  involved  in  a  conical  flow. 

One  can  think  in  terms  of  the  perturbation  being  confined  to 
such  a  small  region  in  the  vicinity  of  the  apex  that  the  dis¬ 
tance  ri  -  r^  itself  is  as  small  as  we  please,  (After  all, 
what  is  a  pointed  cone  from  a  physical  point  of  view,  if  not 
a  blunted  cone  whose  blunt,  irregular,  discontinuous  nose  is 
detectable  only  with  a  microscope?) 

As  a  matter  of  fact,  the  building  up  of  a  conical  flov/ 
arvound  a  blunted  cone,  such  a  flow  being  similar  to  what  should 
be  expected  in  a  pointed  cone,  has  been  investigated  experi¬ 
mentally  (Ref.  12)  and  numerically  (Ref.  13) .  The  conclusion 
was  that  after  about  ten  nose  radii  from  the  tip  of  the  cone  a 
conical  flow  is  well  established.  However,  the  high  entropy 
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layer  consequent  to  the  central  portion  of  the  bow  shock  fills 
a  significant  portion  of  the  shock  layer.  The  entropy  pattern^ 
then,  is  not  what  is  expected  from  a  theoretical  analysis  of 
a  pointed  cone,  and  an  independent  numerical  analysis  is  still 
considered  necessary. 


<lillit.,i. 
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III.  DETAILS  OF  THE  ANALYSIS 

The  Initial  conditions  ■ 

The  present  analysis  differs  from  the  one  mentioned  in 
the  previous  section  (Ref.  13)  in  the  choice  of  the  initial 
conditions  and  in  the  technique  used  to  solve  the  three- 
dimensional  equations  of  a  steady  flow. 

As  said  before,  it  is  convenient  to  choose  the  initial 
conditions  as  close  as  possible  to  the  solution  for  the 
conical  flow,  since  the  computational  time  will  be  reduced 
to  a  minimum.  Moreover,  we  do  not  want  to  start  with  a  set 
of  values  consistent  with  a  blunt  body  calculation,  since  we 
want  to  avoid  the  after  effects  of  a  high  entropy  layer. 

The  initial  values  are  guessed  as  follows: 

It  is  assumed  that  the  trace  of  the  shock  on  a  plane, 
normal  to  the  axis  of  the  cone,  is  an  ellipse  if  the  angle 
of  attack  oc  0.  We  consider  such  a  plane  as  tangent  to  the 
unit  sphere,  whose  center  is  the  apex  of  the  cone,  and  also 
assume  that,  within  the  shock  layer,  0  —  tan  6.  Therefore 
in  Fig.  1,  where  the  plane  is  represented,  the  spherical 
0-coordinate  of  any  point  is  equal  to  the  distance  between 
such  a  point  and  the  origin,  0.  In  this  figure,  O  is  the 
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trace  of  the  axis  of  the  cone  (0=0),  V  is  the  trace  of  the 
vector,  parallel  to  through  the  apex  of  the  cone,  E  is 
the  center  of  the  shock  ellipse,  and  63  are  the  angles 
between  shock  and  body  in  the  leeward  and  windward  side 
respectively,  A  and  B  are  the  semi-axes  of  the  ellipse,  e  is 
the  angle  between  shock  and  body  on  the  x-axis,  and  b  is  the 
semi-angle  of  the  cone.  Note  that 

e  =  A  -  b  (2) 

First,  the  conical  flow  consistent  with  the  assigned 
cone  at  no  incidence  is  computed,  and  a  value,  6^,  of  the 
angle  between  shock  and  body  in  the  plane  of  Fig.  1  is 
obtained.  It  is  then  postulated  that 

(e  +  6a)/2  =  6^  (3) 


independent  of  the  Mach  number  and  the  angle  of  attack  (such 
an  assumption  is  well  justified  by  experimental  results) . 
From  (2)  and  (3)  A  may  be  expressed  as  a  function  of  B  and  C, 
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Consequently,  only  two  parameters,  B  and  C,  are  left  to 
define  the  ellipse.  They  are  related  to  6i  and  63  through 
the  equations, 


B  —  b  +  (6j^  +  63  )/2 
C  =  (6i  -  6a)/2 


(5) 


On  the  plane  of  symmetry,  w  =  0  and  the  derivatives  of  u, 

V,  and  p  with  respect  to  <p  also  vanish.  If  a  conical  motion 
is  assumed,  the  equations  of  motion  in  such  a  plane  reduce 
to  the  simple  form. 


Uq  =  V 


Vg  =  -  u  +  (u  +  V  cot  6  +  w,^  /  sin  0)/(v‘^/a^  -  1) 


v/= 


(6) 


where  a  subindex  means  differentiation  with  respect  to  that 
variable. 

Note  that  the  only  difference  between  (1)  and  (6)  is  the 
presence  of  the  term  w^y4in  0  in  the  latter.  Therefore,  the 
same  technique  used  to  determine  the  conical  flow  at  no 
incidence  could  be  applied  to  solve  (6) ,  if  sin  6  were 


known  as  a  function  of  6.  Since  we  do  not  know  it,  we  assume 


asr* 
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that  w  is  constant  between  M  and  N  and  between  P  and  Q. 

<P 

Let  s  be  the  value  of  6  at  the  shock.  From  the  equa¬ 
tion  of  the  ellipse,  and  letting  9  =  tan  0,  =  1,  it 

follows  that 


s  = 


_  C  cos  <p  +  B  y(B^-C^)/A^  sin^  <P  +  cos^  ^ 


(B/A)  *  sin®  (f)  +  cos®  (p 


and 


^(pcp  =  ■  ^  2B)/a' 


)/a 


(8) 


at  the  symmetry  plane,  where  the  upper  sign  holds  between  M 
and  N  and  the  lov/er  one  between  P  and  Q.  From  the  Rankine- 
Hugoniot  equations. 


(9) 


where  U  is  the  component  of  the  velocity  behind  the  shock 
along  the  normal  to  the  shock  wave,  pointing  inwards,  and 
is  the  derivative  with  respect  to  (p  of  the  component  of  the 
velocity  along  the  tangent  to  the  shock  wave  in  the  plane  of 
Fig.  1.  The  component  of  the  velocity  in  front  of  the  shock 
along  the  normal  is 


V  sin  (s  T  a) 

CO  '  ' 


(10) 
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where  is  the  velocity  at  infinity  and  with  the  upper  sign 
for  point  N  and  the  lower  one  for  point  Q.  From  (ID) ,  U  is 
obtained  through  the  Rankine-Hugoniot  equations.  From  the 
same  equations , 

V  =  V  Cs^^  cos  a  ±  (1  -  s  cot  s)]  (11) 

We  have  now  the  equations  which  ‘are  necessary  and  sufficient 
to  determine  B  and  c.  From  a  computational  point  of  ^^iew, 
an  iteration  procedure  is  convenient  and  is  used.  Two 
values  of  B  and  c  are  guessed,  the  corresponding  values  of 
w  on  the  leeward  and  windward  sides  are  evaluated,  then  the 
system  of  Eqs.  (6)  is  integrated  between  N  and  M  and  between 
Q  and  P.  Two  new  values  of  6i  and  63  result,  from  which 
two  new  values  of  B  and  c  are  obtained,  and  the  procedure  is 
repeated  until  (after  three  or  four  iterations)  the  solution 
is  obtained. 

At  this  stage,  a  number  of  equally  spaced  meridional 
planes  ((p  =  constant)  is  considered  and,  along  the  intersection 
of  each  plane  with  the  unit  sphere,  a  number  of  equally  spaced 
points  is  taken  between  the  cone  and  the  shock,  (left-hand 


side  of  Fig.  2) . 


.4. 
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The  auxiliary  variable  C  is  introduced,  defined  by 


c  = 


8  -  b 


(12) 


where 


6  =  s  ~  b 


(13) 


The  shock  layer  is  thus  mapped  onto  a  rectangle  in  a 
(Y, C) -plane,  where  Y  =  <P-  (right-hand  side  of  Fig.  2).  The 
initial  values  of  u,  v,  w^,  p,  and  p  are  known  at  the  points 
on  the  lines  Y  =  0  and  Y  =  from  the  computation  outlined 
above.  The  shock  shape,  s  =  s{<p)  is  known  as  an  ellipse  and, 
since  we  assumed  that  =  0  (conical  flow) ,  its  geometry  in 

space  is  completely  definqd.  The  Rankine-Hugoniot  equations 
can  be  applied  at  any  point  on  the  shock  (C  =  1) /  to  determine 
the  physical  properties  behind  it,  that  is,  u,  v,  w,  p,  and 
p.  On  the  body  (C  =  0) ,  v  =  0, 

2R  =  R(0,0)  +  R(7r,0)  +  Cr(0,0)  -  R(77,0)]  cos  <p 


2P  =  P(0,0)  +  P(’7,0)  +  Cp(0,0)  -  P(TT,0)]  cos  (p 


(14) 


2w  =  {w^(0,0)  -  +  Cw^(0,0)  +  v/^(7r,0)]  cos  <p}  sin  (p 


where 


R(Y,C)  =  •(.n(p/p„) 
P(Y,C)  =  te(p/p„) 


(15) 


and  u  is  obtained  from  the  definition  of  the  total  enthalpy. 
Note  that,  since  the  entropy,  S,  is  given  by 

S  =  P  -  yR  (16) 

(within  a  constaint  factor,  c^) /  it  follows  that  the  initial 
entropy  on  the  body  is  distributed  with  a  simple  trigono¬ 
metric  law. 

The  values  of  v,  w,  P,  and  R  at  the  interior  points 
of  the  grid  in  the  (Y,C) -plane  are  linearly  interpolated 
between  body  and  shock  and,  again,  u  is  obtained  from  the 
definition  of  the  total  enthalpy. 

At  this  stage,  a  slight  correction  on  the  values  of  p 
(and,  consequently,  of  R)  at  the  interior  points  is  made, 
in  order  to  match  the  total  mass-flow  between  shock  and  body 
with  the  mass-flow  of  the  uniform  stream  through  the  elliptic 
section  of  the  unit  sphere , 
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Let  and  be  the  values  of  p  and  u  as  computed  above,  and 

P  =  pQ  +  C(1  -  O  (18) 

where  T)  is  a  constant,  7)  «  1.  Then  the  mass-flow  between  shock 
and  body  is 

jr  1 

m  =  J  6  ((p)  d<p  I  [pQ  +  1)  ^{1  -  C)  3  (Uq  +  V  )  sin  0  d8  (19) 
o  ‘'o 

where  Ui  is  the  first  order  correction  to  u  as  obtained  from  the 
expression  of  the  total  enthalpy. 


(22) 
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IV.  THE  COMPUTATION  OF  THE  THREE-DIMENSIONAL 

FLOW  FIELD.  INNER  POINTS 

The  computation  of  the  non~conical  flow  subsequent  to  the 
initial  conditions  is  performed  along  the  same  general  lines 
described  for  the  blunt-body  problem  (Ref.  11) .  An  auxiliary 
space  (T,Y,C)  is  used,  where  Y  and  C  have  been  defined  above, 
and  T  =  r.  Referring  to  the  right-hand  side  of  Fig.  2,  one 
can  imagine  a  T-axis  normal  to  the  plane  of  the  figure.  The 
grid  does  not  change  with  T. 

Letting 

q  =  Iniu/Vj 

cr  =  v/u  (23) 

T  =  w/u 

the  equations  of  motion  in  the  new  coordinates  are 

Rt  ^T  ""  "  '''  \ 

+  (Wy  -  Zs^w^/6)/r\i  sin  0  +  (2+0"  cos  0)/r] 

> 

(a*/u®)R^  +  =  -  [Aq^  +  Bq^  -  a*Cs^R^/u®6  -  (a*/yu* )  [  (or/r 

-  Cs^B)S^/6  +  BS^]  -  (V^  +  T^)/2] 

(24) 
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=  -  [Av  +  Bv  +  a®P  /yrufi  +  (v  -  wT  cot  6)/rl 
X  C  Y  ^ 


(a^/yru  sin  0)  (P  -  P^/6) 

i  C  ^  ^  <P  C 


+  (w  +  VT  cot  6) /r] 


[AS^  +  BS^l 


)  (24) 
Cont 'd 


/ 


where 


B  =  r/r  sin  0 


A  =  (a/r  -  Cs^  -  Cs^B)/6 


(25) 


All  points  inside  the  grid  (that  is,  neither  at  ^  =  0  nor  at 
^  =  1)  are  computed  assuming  that 


g(T  +  AT,Y,C)  =  g(T,Y,C)  +  g^(T,Y,C)AT  +  g^,^  (T,  Y,  C)  A  tV2 

(26) 

where  g  is  any  of  the  parameters  g,  v,  w,  R,  and  S.  The  system 
of  Eqs.  (24)  provides  the  formal  expressions  fror  ths  first  derivatives, 
g^.  The  first  derivatives  with  respect  to  Y  and  which  appear 
in  the  right-hand  sides,  are  evaluated  by  central  differences. 

In  addition,  s  is  obtained  from  the  shape  of  the  shock  at  T 
(how  to  obtain  s^  is  explained  in  the  next  section) . 

In  order  to  obtain  the  second  order  derivatives,  g^^,  each 
of  Eqs.  (24)  is  formally  differentiated  with  respect  to  T.  The 
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second  order  mixed  derivatives,  of  the  type  which  appear 

in  the  right-hand  sides,  are  formally  obtained  again  by  differ¬ 
entiating  the  Eqs.  (24)  with  respect  to  ^  and  Y.  In  this  way, 
the  evaluation  of  the  last  term  in  (26)  is  reduced  to  computing 
3-point  centered  second  derivatives  with  respect  to  ^Y,  and 
YY,  and  to  some  algebraic  manipulations.  Note  that  6,j,  = 

The  terms,  s  s  „  are  evaluated  eis 
rT  (pT 


s^(T  +  A  T,Y)  -  s^(T,Y) 


s  (T  +  A  T,Y)  -  s  (T,Y) 

_ 52 _ 

'  A  T 


Derivatives  ivith  respect  to  ^  and  Y  of  complicated  factors,  such 
as  A,  a^/yru  sin  0,  etc.  are  not  formally  evaluated,  but  com¬ 
puted  by  central  differencing  of  the  factors  themselves. 

The  value  of  A  T  which  can  be  used  at  each  step  is  limited 
by  a  stability  criterion.  The  well  known  Couiant-Friedrichs-Lewy 
rule  is  applied,  which  forbids  computing  values  at  a  point  outside 
the  domain  of  influence  of  the  initial  points  (Ref.  16) . 
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V.  COMPUTATION  OF  THE  SHOCK  WAVE 

The  points  on  the  shock  wave  are  computed  with  an  improved 
version  of  a  technique  successfully  used  in  general  three- 
dimensional  flow  field  computations  (Ref.  13)  . 

The  equations  of  motion  in  spherical  coordinates  are  written 
moving  all  derivatives  with  respect  to  (p  to  their  right-hand 
sides,  and  the  third  momentum  equation  is  stricken  out,  so  that 
the  system  looks  as  follows: 

-  S^/y  -  aSg/yr  +  P^/y  +  aPg/yr  +  q^  + 

+  w  /u)/sin  6  +  2  +  a  cot  9]/r 


(a^/u^)  (-S^y  -  crSg/yr  +  P^/y)  +  +  crq^/r  =  [  (a^S^/yu® 

-  q  )  T/sin  6  +  CT®+  T^]/r 

a^P_/yru  +  v  +  crv^/r  =  [-v  -  Tv  /sin  0  +  wt  cos  61/r 
0  r  0  (p 

S  t  CTSfl/r  =  -  TS  /r  sin  6 
r  0  <p 


(27) 


In  a  (r, 6) -plane,  this  system  has  three  characteristics,  defined 
by 


\ 


I 


d8  _  uv  +  a  /y/u^+v^-a^ 
dr  ~  u^  -a® 


^  II  _  V.  _  uv  -  a  /y/u^+v^-a^ 
~  dr  ”  u^-a® 


dr 


cr 


(28) 
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The  compatibility  equation  along  the  first  characteristic  reads 


(a/yu=  )  +v=  -a=  +  ^  = 

dr  dr 


r 


(w  /u  sin  0+2 
<P 


+  or  cot  6)  +  X  (rq  /sin  0-cr^-  T® ) 

<P 


-  (rv^/sin  0  +  V  -  wT  cos  0)/ru 

-  [aq  +  v^^+v^-a2  |  p  ra/yr  (u® -a® )  sin  0 

(29) 

By  computing  the  coefficients  and  the  right-hand  side  at  a 
suitable  initial  point,  the  values  of  P  and  q  -.ay  match  the 
values  of  P  and  q  computed  from  the  Rankine-Hugoniot  conditions, 
with  a  proper  choice  of  s^.  By  iteration,  the  matching  values, 
including  s^,  are  easily  computed. 


VI .  COMPUTATION  OF  POINTS  ON  THE  BODY 

A  similar  technique  is  used  to  compute  points  on  the  body. 
The  compatibility  equation  along  the  second  characteristic, 


(a/yu®) 


(w  /u  sin  6+2 


+  ff  cot  0)  “  (Tq^/sin  0-y^ -T® ) 


+  (tv  /sin  0  +  V  -  wT  cos  0)/ru  + 

(fi 


P  ra/yr  (u^ -a^  )  sin  0 
‘P 

(30) 


must  be  used.  Again,  the  coefficients  and  the  right-hand  side 


are  computed  at  a  suitable  initial  point  and  the  condition, 

V  =  0  on  the  body,  is  used.  Consequently,  P  is  found  at  a  body 
point.  The  fluid  momentum  equation  then  is  used,  in  which  v  =  0 


1  Dw  T 

~  —  =  w  +  -  w 

w  Dt  r  r  sin  0  0 


(w  +  a®P  /yu  sin  6)/r  (31) 


the  last  of  Eqs.  (27)  reads  now 


1  ^ 
u  Dt 


+ 


T 

r  sin  0 


S 

<P 


0 


(32) 


Both  Lagrangian  derivatives  are  taken  along  a  streamline  which 
wets  the  body.  Thus,  by  finding  the  initial  point  on  such  a 


streamline  and  computing  the  pertinent  values  of  w,  S,  and  the 


right-hand  side  of  (31) ,  the  values  of  w  and  S  at  the  point  on 
the  body  may  be  evaluated.  Finally,  the  value  of  u  at  the  body 
point  is  obtained  from  the  condition  of  constant  total  enthalpy. 

At  the  body  point  in  the  leeward  symmetry  plane  (<p  =  0)  , 
however,  the  entropy  is  forced  to  be  equal  to  the  entropy  of 
the  neighboring  point  in  the  same  plane.  A  similar  policy  is 
adopted  for  the  body  point  in  the  windward  symmetry  plane 
((/)  =  it)  .  The  reason  for  it  is  that  no  streamline  enters  or 
leaves  the  symmetry  planes.  When  the  flow  has  settled  down  to 
a  conical  one,  the  entropy  on  each  of  these  planes  must  be  a 
constant.  Our  device  allows  the  body  points  in  the  symmetry 
planes  to  adjust  their  entropies  to  the  final  entropies  at  the 
shock  points,  independent  of  the  original  guess. 
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VII.  RESULTS  AND  DISCUSSION 

The  following  sample  problems  have  been  run,  in  order  to 
test  the  technique: 


(1) 

M  =  7.95,  €  =  10°, 

CO 

a  = 

8° 

(2) 

M  =  7.95,  C  =  10°, 

CO 

a  = 

4° 

(3) 

00 

II 

8 

S 

€  =  20°, 

a  = 

5° 

(4) 

M  =  20, 

00 

C  =  30°, 

a  = 

10° 

The  first  two  cases  were  chosen  to  compare  our  results  with 
Tracy's  experiments  (Ref.  14),  the  third  one  to  conclude  the  dis¬ 
cussion  and  to  check  with  other  experiments  which  were  made  in 
testing  our  three-dimensional  flow  field  computations  (Ref.  13) , 
the  fourth  one  to  show  the  reversed  situation  (leeward  side  thinner 
than  the  windward  side)  which  is  to  be  found  at  high  Mach  numbers 
(Ref.  15) .  We  will  analyze  the  first  case  in  detail. 

It  took  four  iterations  to  define  the  shape  of  the  initial 
shock  wave.  Table  I  shows  the  values  of  B  and  C,  their  errors 
(that  is,  the  differences  between  B  and  C  at  one  step  and  B  and  C 
at  the  next  step)  and  the  corresponding  values  of  w^  in  the 
leeward  side  (labeled  1)  and  the  windward  side  (labeled  2) .  In 
the  same  table,  values  of  u,  v,  p,  p,  M,  and  v/u  are  given  as 
functions  of  0  in  the  symmetry  plane  (from  the  shock  to  the  body) . 
Note  that  all  dimensional  quantities  are  made  nondimensional; 


23 


pressures  and  densities  are  divided  by  and  respectively; 
velocity  components  are  divided  by  •  The  subsequent 

computation  was  performed  (in  this  as  in  the  other  cases)  con¬ 
sidering  11  meridional  planes,  spaced  18°  apart,  and  6  points 
in  each  plane  (including  the  shock  and  body  points) . 

As  a  result  of  the  non-conical  computation,  the  shock 
v/ave  undergoes  a  deformation.  The  values  of  s^,  which  are 
assumed  as  zero  initially  and  should  be  identically  zero  in  a 
conical  flow,  change  with  r  and  (p.  If  the  flow  sets  down  to 
a  conical  one,  uhen  s^  must  vanish  again  identically.  A  first 
indication  of  the  tendency  of  the  flow  toward  a  conical  state 
may  be  found  by  plotting  the  range  of  values  of  s^  at  each  com¬ 
puted  value  of  r.  This  is  done  in  Fig.  3.  Since  s^  is  of  the 
order  of  .2,  the  shock  shape  may  be  considered  as  stabilized 

when  s  is  of  the  order  of  10""^ ,  and  this  situation  occurs  at 

r 

step  200,  r  «  13.6,  far  before  the  end  of  the  run,  which  was 
interrupted  at  step  540,  r  =  1073. 

In  Fig.  4  the  internal  mass-flow  (across  a  spherical  sec¬ 
tion  of  the  shock  layer  at  a  given  station)  and  the  external 
mass-flow  (across  the  whole  shock  wave  from  the  apex  to  the  same 
station)  are  plotted  after  having  been  divided  by  r® .  It  is 
clear  again  that  a  stable  flow  configuration  is  reached,  after 
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about  400  steps  (in  a  conical  flow,  m/r®  is  a  constant) .  It  is 
also  interesting  to  note  that  the  error  (difference  between 
internal  and  external  mass  flow)  is  of  the  order  of  1/1000  of 
the  total  mass  flow) . 

There  is  a  way  of  testing  the  soundness  of  the  results. 

If  the  flow,  as  found  by  the  present  technique,  is  conical,  then 
it  is  ruled  by  Eqs.  (6)  in- the  symmetry  plane.  Let  us  perform 
again  the  integration  of  such  a  system,  starting  from  the  shock 
as  defined  by  the  present  technique  and  using  what  we  assume 
are  the  correct  values  of  ^^(9) •  The  integration  should  give 
the  correct  value  for  the  body  angle  and  the  distributions  of  u 
and  V  between  shock  and  body  should  match  the  ones  obtained  with 
the  present  technique.  This  actually  happens  within  three  sig¬ 
nificant  figures.  The  body  angle,  in  radians,  is  .17453.  The 
body  angles  computed  from  (6)  with  the  values  of  w  interpolated 
from  the  last  step  in  our  sample  run  are  .17410  in  the  leeward 
side  and  .17445  in  the  windward  side.  The  distribution  of  v(0) 
as  computed  with  the  present  technique  and  from  the  integration 
of  (6)  is  shown  in  Fig.  5. 

A  word  about  the  entropy  distribution  is  in  order  here.  As 
pointed  out  by  Ferri  (Ref.  7),  the  constant  entropy  lines  on  the 
unit  sphere  converge  to  a  singular  point  on  the  leeward  plane  of 
symmetry.  In  our  computation,  the  initial  (arbitrary)  distribution 
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of  entropy  is  the  one  shown  in  Fig.  6,  left-hand  side.  At  the  end 
of  the  computation,  the  constant  entropy  lines  appear  as  in  the 
right-hand  side  of  the  same  figure.  The  singularity  has  been 
built  up  throvigh  a  numerical  procedure.  Note  that  such 'a  result 
is  obtained  using  a  very  coarse  mesh.  To  stress  the  point, 

Fig.  7  shows  the  building  of  the  singularity  by  plotting  several 
distributions  of  entropy  on  the  body,  at  different  steps  in  the 
computation. 

The  pressure  distribution  is  shown  in  Fig.  8.  In  Fig.  9  the 
pressure  distribution  on  the  body  is  plotted  as  a  function  of  (f>  for 
case  (3)  specified  above.  As  expected,  no  sizeable  difference 
exists  among  the  present  results,  those  obtained  by  a  proper  usage 
of  Kopal's  tables  (Ref.  2),  and  the  pressure  distribution  on  a 
blunted  cone  at  a  certain  distance  from  the  nose.  The  latter 
results  have  been  obtained  analytically  (by  the  technique  of  Ref.  13) 
and  experimentally  (Ref.  12). 


~ 
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VIII.  LIMITATIONS  OF  THE  INVISCID  MODEL 

To  conclude,  Pig.  10  shows  the  shock  shapes  for  the  four 
cases  computed.  In  the  first  two  cases,  some  points  on  the 
shocks  as  determined  experimentally  (Ref.  14)  are  shovm  by  dots. 
The  agreement  between  computation  and  experiment  in  the  windward 
side  is  a  proof  of  the  soundness  of  the  numerical  procedure.  The 
disagreement  in  the  leeward  side  was  expected.  The  present  model 
is  one  of  an  inviscid  flow.  In  reality,  the  viscous  effects  are 
important  on  a  cone  at  an  angle  of  attack,  due  to  the  three- 
dimensional  behavior  of  the  boundary  layer  (Ref.  17) .  By  com¬ 
parison  of  the  numerical  results  with  other  data  from  Ref.  14, 
it  appears  that  the  difference  in  thickness  in  the  leeward  sym- 
metry  plane  is  of  the  order  of  the  thickness  of  the  viscous  core. 

In  conclusion,  the  present  computation,  although  self- 
consistent  and  accurate,  is  not  yet  adequate  to  describe  the 
real  motion,  for  lack  of  a  realistic  model  of  the  gas.  A  more 
complicated  model,  in  which  viscous  effects  are  taken  into 
account,  will  be  presented  in  a ^ further  communication. 
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TABLE  I 


COXE  AT  AN  ANGLE  OF  ATTACK 
PROGRAH  3A  -  RUN  NUNBER  It  ON  ll/li/65 


NA«  5,  MA«  IC.  KA-  2jr,  JA»  10.  LA«  -0.  IB«  -0,  LC«  -0,  LO*  -0.  LE«  -C,  LF«  -0 

I 

CGNE  HALF  ANGLE*  1.03399F  01  DEGREES.  ALPHA*  S.OOOTOE  00  DECREES.  UO*  9.410C0E  00.  CANHA*  l.ACOOOE  00.  STAB*  5.C''CC0E  r 
GUESS  FOR  SHOCK*  1.5C0P0F  01.  0ELTH*lC.C0CCCE-e3 
EPS  I  TO  7 

10.0300E-05  lC.OCOOE-05  lO.OOOOE-05  l.OOOOE  00  -0.  -0.  -0. 


NACH  NUNRER*  7.9529E  00.  STANDOFF  DISTANCE  AT  2ERQ  ANGLE  OF  ATTACK*  5.297662E-02 


SHOCK  ITERATION 


B 

C 

ERRtai 

ERR(C) 

riPHtll) 

MPHl 12) 

1 

2.461303E-01 

3.321745E-02 

1.86236aE-02 

3.321 f45E-02 

-1.309597E 

00 

1.309597E  00 

1 

2.495162E-CI 

2.7(IS626E-02 

3.385855E-C3 

-5.351U  >03 

-1.253872E 

Of 

6.1T2868E-C1 

2 

2.494329E-01 

2.8029SAE-02 

-8.35556SE-04 

-5.79149CE-04 

-1.265693E 

00 

6.350081E-01 

J 

2.493270E-ei 

2.818341E-C2 

-1.223702E-04 

1.340313E-04 

-1.26673SE 

00 

6.52B740E-01 

4 

2.493209E-01 

2.813864E-02 

1.199171E-05 

3.414042E-05 

-1.266764E 

00 

6.527720E-01 

GUESS  FOR  LEEMARO 

SHOCK 

ENTROPY* 

3.507286E-04 

THETA 

U 

V 

RHO 

P 

v/u 

1 

2.77S095E-31 

9.320688E 

00 

-1.117585E  OC 

1.157315E 

00 

.1,227396E 

03 

7.T94015E 

00 

-1.199037E-01 

2 

2.675095E-01 

9.331864E 

00 

-9.95993JE-01 

1.176324E 

DO 

1.255712E 

00 

7.6  76840E 

Of 

-1.067304E-01 

3 

2.575C95E-01 

9,341824E 

00 

-8.855479E-01 

1.184930L' 

00 

1.268592E 

00 

7.664T06E 

00 

-9.479391E-02 

4 

2.475O95E-01 

9.3506r9E 

00 

-T.781674E-01 

1.190135E 

90 

1.276401E 

00 

7.65r4l8E 

OC 

-8.321828E-02 

5 

2.375095E-01 

9.3S84S0E 

00 

-6.720345E-01 

1.193427E 

00 

1.281347E 

00 

7.6528288 

OC 

-7,18103TE-02 

6 

2.275r95E-0l 

9.365181G 

00 

-5.664931E>01 

t.l95387E 

00 

1,284294E 

00 

7.650102E 

00 

-6.048469E-02 

7 

2.175095E-01 

9.3708A5E 

00 

-4.609509E-01 

1.196286E 

00 

1.285646E 

00 

7.648854E 

00 

-4.918993E-02 

8 

2.075095E-01 

9.375455E 

00 

-3.552157E-01 

1.19624tE 

00 

1.285578E 

00 

7.648916E 

30 

-3.78878AE-02 

9 

1.975095E-01 

9.379007E 

00 

“2.A89646E-01 

1.195278E 

00 

1.284130E 

CO 

7.650253E 

00 

-2.654488E-02 

10 

1.8T5O95E-01 

9.381696E 

00 

-l.419U3e-0l 

1.19334SE 

00 

1.281227E 

00 

7.652938E 

00 

-1.512705E-02 

11 

l,775C95E-0l 

9.382915E 

00 

-3.3T4197E-02 

1.190322E 

00 

1.276682E 

00 

7.657156E 

00 

-3.59610TE-03 

12 

1.7S0C95E-C1 

9.333300E 

00 

-5.4891T7E-03 

l.)90135E 

00 

1.276400E 

00 

7.557418E 

30 

-5.85C13CE-04 

13 

1.74697OE-01 

9.3830CIE 

30 

-1. 9120788-03 

1.190133E 

00 

1.276397E 

00 

7.657421E 

OP 

-2.037810E-04 

14 

1.745407E-01 

9.383032E 

00 

-1.213769E-04 

t.l9')132E 

00 

1.2T6396E 

OP 

7.657422E 

90 

-1.2935836-05 

15 

1.745S10E-O1 

9.3B3002E 

00 

-9.362918E-C6 

1.190132E 

00 

i.2T6396E 

OC 

7.657422E 

00 
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